Dynamical probing of a topological phase of bosons in one dimension 
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We study the linear response to time-dependent probes of a symmetry-protected topological phase 
of bosons in one-dimension, the Haldane insulator (HI). This phase is separated from the ordinary 
Mott insulator (MI) and density-wave (DW) phases by continuous transitions, whose field theoretical 
description is here reviewed. Using this technique, we compute the absorption spectrum to two 
types of periodic perturbations and relate the findings to the nature of the critical excitations at 
the transition between the different phases. The HI-MI phase transition is topological and the 
critical excitations possess trivial quantum numbers: they correspond to particles and holes at zero 
momentum. Our findings are corroborated by a non-local mean-field approach, which allows us to 
directly relate the predicted spectrum to the known microscopic theory. 



I. INTRODUCTION 

According to the Landau paradigm, any two phases 
with the same symmetries can be adiabatically connected 
through some (complicated) path in parameter space (the 
best known example being the liquid and gas phases of 
water). For many years, the only counter-example to 
this paradigm was the quantum Hall states, which do 
not break any local symmetry, but nevertheless cannot 
be connected together without crossing a phase transi- 
tion. In recent years, this phenomenon has been extended 
to a broader class of phases, known as "topological band 
insulators'*^. Lacking a local order parameter, topolog- 
ical insulators may possibly be used to realize quantum 
circuits unaffected by local sources of dissipation, such 
as disorder and time-dependent noised However, for the 
same reason, their experimental characterization is gen- 
erally very hard: local probes are unable to distinguish 
between the bulks of a conventional insulator and of a 
topological insulator. 

This problem is usually overcome by studying inter- 
faces between topological phases and trivial phases (such 
as the vacuum), where protected edge states must ap- 
pear. In the case of the integer quantum Hall states, 
the number of edge modes coincides with the number of 
filled Landau levels and leads to a quantized transverse 
conductance. Similarly, the appearance of topologically- 
protected zero-energy modes at the boundaries between 
different phases constituted the first experimental evi- 
dence of the topological insulatoripHSl A second strat- 
egy consists in measuring non-local "string" order pa- 
rameters, characteristic for example of topological phases 
in one dimension. A proof-of-principle experiment has 
already been successfully demonstrated using ult racold 
atoms in optical lattices^, following our prediction d 10 ' 11 '. 

In the context of ultracold atoms, a third alternative 
method consists in probing the dynamic response of the 
syste m in t he vicinity of a transition between different 
phasea^ES. Here, we specifically refer to continuous 
phase transitions, which can be detected by the closing of 
a gap in the response to a local probe^H In the vicinity 
of these transitions, the low-frequency excitations hold 
information about the nature of the transition. In par- 



ticular, because transitions between topologically trivial 
and non-trivial phases do not break any local symme- 
try, their low-frequency excitations must possess trivial 
quantum numbers. In contrast, if a transition breaks a 
particular symmetry, its low-frequency excitations must 
possess the same quantum number as the order param- 
eter of the symmetry broken phase. Dynamical probes 
couple to these low-frequency excitations and may pro- 
vide useful information about the nature of the transi- 
tion. 

In the case of translational invariance, for example, at 
the transition from a homogeneous phase to a density- 
wave phase, the critical excitation must possess the same 
momentum as the density-wave order parameter. In con- 
trast, at a topological phase transition the critical exci- 
tations always possess zero momentum. Similarly, sys- 
tems that conserve the total number of particles possess 
an additional U{1) "gauge" symmetry, which is sponta- 
neously broken at the transition to a superfluid phase. 
At this transition the critical excitations are Bogoliubov 
quasi-particles, which do not conserve the total number 
of particles (or to be more precise, shift the number of 
particles by a non-integer amount). On the contrary, at 
a transition between a normal and topological phase, the 
symmetry is preserved and, accordingly, the critical ex- 
citations must carry integer charge. 

In this article we consider a particular example of 
one dimensional symmetry- protected topological insula- 
tor, the "Haldane Insulator' ™ ' 11 ' 15 ' 16 !. This phase arises 
in models of interacting bosons on a one-dinrensional lat- 
tice, which are relevant to ultracold dipolar atoms^ or 
molecules^ and Rydberg atomgED, in optical lattices^. 
These systems support additional phases that are topo- 
logically trivial, such as the Mott insulator (MI) and the 
density-wave (DW) phase. In the simplest DW phase 
the density oscillates between even and odd sites, spon- 
taneously breaking the discrete translational invariance, 
modulo two lattice sites. As we will see, the dynamic 
response demonstrates that the critical excitations of the 
HI-DW transition are phonons at momentum k = n. In 
contrast, at the topological MI- HI transition the dynamic 
response demonstrates that the critical excitations are 
simply particles and holes at zero momentum. 



This article is organized as follows. In Sec. [TT]we de- 
scribe the ground state properties of the Haldane insu- 
lator and present new numerical results aimed to char- 
acterize its low-lying excitations. In Sec. |III| and [TV] we 
introduce the field theoretical description of the Haldane 
insulator and use this method to predict the dynamic re- 
sponse of the system. Our findings are in contrast with 
the predictions of Ref. [TO1 obtained through a non-local 
mean field approach. In Sec. [V] we show how to cor- 
rectly apply the mean-field approach to obtain results 
consistent with the field-theoretical description. Sec. |VI| 
concludes the article with a summary of the results and 
an outlook to future research directions. 



II. THE HALDANE INSULATOR 

Lattice bosons with non-local interactions can form 
a zoo of different phases with local order parameters, 
including superfluids (SF), density- waves (DW), super- 
solids, and dimerized phases (See for example Ref.s l2lT 
25]) . Traditionally, the only phases without a local pa- 
rameter were Mott insulators at integer filling. These 
phases are topologically trivial because their ground 
states can be adiabatically deformed to site-factorizable 
states with an exact (integer) number of atoms per site 
( |11111...), |22222...), ... states ) In this sense, the Mott 
insulators (MI) of bosons are analogous to the trivial 
band insulators of fermions. In Ref. [HJ1 we found that 
non-local (short-ranged) interactions can stabilize a non- 
trivial topological phase, the Haldane insulator (HI). 

In our theoretical study, we considered the one- 
dimensional Bose-Hubbard model with extended inter- 
actions: 
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The first term t describes the tunneling ( "hopping" ) be- 
tween neighboring sites, the second term U describes the 
effective interactions between atoms on the same site, 
and Vj describes non-local interactions. This model ap- 
plies to ato ms and molecules with fixed electric or mag- 
netic dipole ^' l 18 * 24 !, and to Rydberg atoms with induced 
dipoles^l. For simplicity, we will assume that only the 
nearest-neighbor term V\ = V is non zero. This assump- 
tion is not essential for the present analysis, provided 
that Vj decays fast enough. 

At V = the Bos e-Hub bard model describes a quan- 
tum phase transitio n 121 * 25 * * 2 ^ between a one-dimensional 
superfluid (for small U/t) and a Mott insulator (for large 
U/t) . In the presence of a finite non-local repulsion 
V > (of the order of U) the Hamiltonian (fl) de- 
scribes a tran sition to a topological phase, the Ha dane 
insulato j 10 ! 11 !. This transition was first discovered by nu- 
meric studies of the excitation gap, reproduced in the 
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FIG. 1. Phase diagram of the extended Bose Hubbard model 
(II]) at filling of fl — 1 particles per site, as obtained by Density 
Matrix Renormalization Group (DMRG) calculations with 
imax = 256 sites and M max = 250 states per block. IN- 
SET: Excitation gaps A c and A n (defined in the text) along 
the path U/t — 6. The phase transitions are located at the 
points in which the gap vanishes (V/t » 4 and V/t « 4.75). 



inset of Fig[T] This graph clearly shows the existence of 
three gapped phases, separated by two continuous phase 
transitions. The three gapped phases are the Mott insu- 
lator (MI), the Haldane insulator (HI), and the density- 
wave (DW) phases. This last phase is adiabatically con- 
nected to a classical state with an alternating number of 
bosons per site (| 202020...)), and spontaneously breaks 
the lattice translational invariance, modulo two sites (Z^ 
symmetry) . 

The Hald ane insu lator is a symmetry-protected topo- 
logical phase^E^EGD j n analogy to the topological insu- 
lators of fermions (which are protected by time-reversal 
symmetry), the HI is protected by an anti-unitary dis- 
crete symmetry, the lattice-inversion symmetry. In the 
presence of this symmetry, the HI cannot be adiabati- 
cally connected to a site-factorizable state. Remarkably, 
the lattice-inversion symmetry is trivially broken at the 
edges of the lattice, and therefore the Haldane insulator 
does not have protected edge states^. 

Because both the Mott insulator and the Haldane in- 
sulator do not break any local symmetry, they cannot 
be distinguished by any local order parameter. Exploit- 
ing an analogy to spin chains, we demonstrated that 
these two phases can nevertheless be characterized by two 
"string" order parameters^*. These order parameters are 
non-local in the sense that they consist o f an in finite prod- 
uct of operators acting on different site d 30 * 31 *. Their ex- 
perimental detection therefore requires the simultaneous 
measurement of all the sites in the lattic e. Usi ng new mi- 
croscopic techniques for ultracold atoms^ES the string 
order of the Mott insulator has been recently measurecP. 



The existence of a string order is however specific to 
one dimension and, to the moment, does not seem to ex- 
tend to other topological phases^ . In this article, we 
present an alternative method to characterize the topo- 
logical nature of these phases. Our approach consists in 
probing the low-energy excitations of the different phases. 
It is therefore not specific to one-dimension and can be 
easily generalized to higher dimensions. In one dimen- 
sion, both methods apply and could be used to comple- 
ment and verify one-another. 



Elementary excitations 

As explained in the introduction, we propose here to 
study the topological properties of the Mott and Haldane 
phases by probing the low-energy excitations at the tran- 
sition between the two phases. To characterize their na- 
ture, let us consider again the inset of FigJTj This graph 
shows two independent curves. One curve (diamonds) 
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represents the "neutral gap" A„ = E 
which is the gap to the first excited state with the same 
number of particles as in the ground state (8N = 0). It 
was calculated by targeting the first excited state in the 
DMRG calculation. The second curve (squares) repre- 
sents the "particle-hole gap" 2A C = E SI ^ =1 + E s J r= _ 1 — 

2Egjy =Q , which was calculated by adding or subtracting 
one particle (SN = ±1) and computing the new ground 
state. 

At the MI-HI phase transition the particle-hole gap 
closes, indicating that at this point the energy needed 
to create a particle or a hole ("the mass term", in the 
field-theoretical terminology) becomes zero. At the HI- 
DW phase transition, on the other hand, the particle-hole 
gap remains open, indicating that the phase transition is 
due to a different type of excitations. We can easily guess 
the nature of these excitations by recalling that the HI- 
DW phase transition involves the spontaneous breaking 
of a local i?2 symmetry, related to translational invariance 
modulo two lattice sites. In order to break this symmetry, 
the ground state has to mix with a phonon at the edge 
of the Brillouin zone (i. e. at momentum k = ir), whose 
excitation gap must be closed at the phase transition. 

To verify this hypothesis, we compute the overlap be- 
tween the lowest excited state in the vicinity of the 
HI-DW phase transition and a trial wavefunction for a 
phonon at momentum k: 
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Here | ^ gs ) is the ground state wavefunction as given 



by the DMRG algorithm and Srij 



the difference 



between the local occupation and the average occupation 
n. As shown in Fig. [5J the numeric result confirms the 
existence of a large overlap with a phonon at momentum 
k — 7r, supporting our assumption. 




FIG. 2. Characterization of the lowest neutral excitation 
*?2_ ) ■ (a) Matrix element of M» = ( *^ =0 | 5m | * gs ) 
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(b) Square of the Fourier transform 

M k = £.exp(tfeCi/L).Mi = (* (1) 

of momentum k. The pronounced peak of \M k 

confirms that | ^jv = o ) nas a ^ ar S e overlap with a phonon at 
the edge of the Brillouin zone 



III. CRITICAL THEORY OF THE PHASE 
TRANSITIONS 

The properties of any system near a continuous phase 
transitions do not dependent on the details of the micro- 
scopic model and are determined by a universal critical 
theory. The critical theories of the MI-HI and HI-DW 
phase transitions were derived in Ref. [TT] using a field- 
theoretical approach, the "bosonization of bosons"'2212SI. 
In what follows, we outline the properties of these theo- 
ries, which will be used to compute the dynamic response 
of the system in Sec. |IV| 



MI-HI phase transition 

The critical theory of the MI-HI phase transition is 
described by the sine-Gordon model 



H = dx 
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Here the field 4>{x) measures the displacement of the par- 
ticles with respect to a fixed lattice and is related to 
the particles' density by n(x) ~ n — (l/ir)d x (f>, where 
uq = n/a, n is the average number of particles per site, 
and a is the lattice constant. The field 0(x) is canonically 
conjugate to n{x) and corresponds to the coherent phase 
of the bosons. 

The quadratic part of the Hamiltonian pi simply de- 
scribes linear phonons, i.e. the Goldstone mode associ- 
ated with translational invariance, in one dimension. The 



unitless prefactor K is non-universal and it varies along 
the MI-HI phase boundary. A renormalization group 
(RG) analysis^ shows that, in order for a direct MI-HI 
phase transition to occur, K must be bounded between 
1/2 and 2. For K > 2, the MI and HI phases would be 
separated by a superfluid phase, and for K < 1/2, by a 
third gapped phase. 

The cosine term in Eq. (|3]) describes the effects of the 
periodic potential (lattice). This term pins the field <j>(x) 
to a fixed value, modulo n. The MI-HI phase transition 
corresponds to a change in the sign of the prefactor. In 
the Mott insulator g > and the field is pinned to zero, 
indicating that the particles are centered at the lattice 
sites. In the Haldane insulator g < and c6 is pinned 
to 7r/2. The particles are centered at the lattice bonds, 
and live in a coherent superposition of two (or more) 
neighboring sites. A weak-coupling analysis leads to g s=s 
Pn(U — V), indicating that the MI-HI phase transition 
should be located approximately at U = V . 

The low-energy excitations of the MI-HI phase transi- 
tion correspond to jumps of ±7r in the field <fi(x), termed 
"kinks" , or "solitons" . Because the density of particles 
is proportional to the derivative of c6, each kink involves 
the addition, or subtraction, of one particle. Using a per- 
turbative RG approach it is possible to show that each 
kink increases the energy of the state by a charge gap 
A c ~ ^VC 2 — -*0. In addition, depending on the value of 
K, the kinks may either repel (for K > 1), or attract (for 
K < 1). In the former case, the neutral excitation gap 
is given by twice the energy cost of a kink, A n = 2A C . 
In the latter case, particles and holes form bound states 
("breathers"). Their total mass is lowered by the bind- 
ing energy of the pair. In what follows we consider only 
1 < K < 2, corresponding to intermediate values of the 
lattice depth. 

The sine-Gordon model ^ is a useful representation 
of the critical theory of the MI-HI phase transition be- 
cause it is exactly solvable. In particular, the matrix 
elements between the ground state and any excited state 
can be obtained following the "form factor" approach 
introduced by Smirnov^S This approach relies on the 
exact expression for the scattering-matrix between kinks 
and breathers. 



HI-DW phase transition 

The transition from the Haldane insulator (HI) to the 
density-wave (DW) phase breaks a local Z^ symmetry. In 
Ref . ITT1 we showed that, indeed, this transition belongs 
to the universality class of the Ising transition in 1+1 
dimensions. The critical theory of this transition is cap- 
tured, for instance, by the spin-1/2 transverse-field Ising 
chain 



Here erf and af are Pauli matrices. The model Q shows 
a quantum phase transition at A = 1 between a para- 
magnet, with zero magnetization (erf) = 0, and a fer- 
romagnet, with (erf) ^= 0. The ferromagnet breaks the 
Zi symmetry, related to the rotation of 7r around the z- 
axis. The transverse- field Ising model (4| has an exact 
solution in terms of fr ee ferm ions, based on the Jordan- 
Wigner transformation^^. This mapping shows that 
the excitation gap vanishes linearly and symmetrically 
around the phase transition: 



A=|l-A| 



(5) 



The exact solution also shows that the nature of the 
low-energy excitations changes dramatically across the 
phase transition. In the paramagnetic phase, the low- 
energy excitations are local spin-flips, created by the 
operator af. In the ferromagnetic phase, on the other 
hand, the elementary excitations are domain walls be- 
tween spin-up and spin-down regions. 

In the original Bose-Hubbard model the paramagnetic 
and ferromagnetic phases correspond respectively to the 
Haldane insulator and the density-wave phase. The spin- 
up and spin-down states correspond to the two possi- 
ble configurations of the density-wave (|202020...) and 
1 020202...)), leading to the mapping 
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Eq. ((6) indicated that the low-energy excitations of 
the Haldane insulator are density fluctuations (phonons) 
with momentum n. In contrast, the low-energy exci- 
tations of the density-wave phase are non-local domain 
walls. 



IV. DYNAMIC PROBES: ABSORPTION 
SPECTRA 



We now apply the above field-theoretical description 
to compute the dynamic response of the system in the 
vicinity of the MI-HI and HI-DW phase transitions. In 
particular, we focus on two specific perturbations: lattice 
modulations and Bragg spectroscopy at momentum k = 
7r. As we will see, these perturbations arc easily realizable 
with ultra-cold atoms and best capture the differences 
between the two transitions. The results obtained by 
this method are summarized in Fig.s [3] and [4] 

In an idealized version of the proposed experiment the 
system is initially prepared in the ground state, a periodic 
perturbation is applied with frequency u) for a given time 
£pcrt 3> 1/w, and finally the amount of absorbed energy 
dE is computed, by measuring the increase in tempera- 
ture of the system. The absorption spectrum is defined 
as the ratio between the absorbed energy and the pertur- 
bation time ipcrt- According to Fermi golden rule, this 



ratio is independent on i per t, and equals to 
S(u) = J2\(^\T\V gs )\ 2 (e a - e gs ) S(lu + e gs - e a ). 

a 

Here | \l/g S ) and | ij) a ) are, respectively, the ground- and 
excited-states of the unperturbed Hamiltonian, with en- 
ergies e gs and e a (we work in units where H = 1). The 
operator T models the perturbation and will be derived 
below for two specific cases, denoted by Th op and Tfc =7r . 

The first perturbation consists of small modulations of 
the lattice strength. It mainly influences the tunneling 
matrix element between two neighboring sites and corre- 
sponds to the operator 
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The second perturbation corresponds to the addition of 
a superlattice with a given momentum k, which locally 
modifies both the tunneling and the chemical potential. 
In the case of k — it, the tunneling is modulated uni- 
formly along the lattice, while the chemical potentials of 
odd and even sites are shifted in opposite directions. The 
correspondent operator is the sum of fl8| and 



T fc=7r = ^y^(-i)'fa, 



(9) 



Bragg spectroscopy, on the other hand, does not couple 
to the low-energy excitations of the MI-HI phase transi- 
tions. The critical excitations of this phase transition 
are particles and holes with zero momentum, while Tfc =7r 
couples only to excitations with total momentum k = tt. 
Thus, Bragg spectroscopy at momentum k — it does not 
show any particular low-frequency feature at the MI-HI 
phase transition, despite the closing of the gap. This ob- 
servation will allow us to clearly distinguish between this 
transition and the HI-DW transition to be considered be- 
low. 



HI-DW phase transition 

The critical theory of the HI-DW phase transition is 
captured by the transverse-field Ising model Q . The op- 
erator Thop corresponds to the term of the Hamiltonian 
(Til) which stabilizes the Haldane insulator, and therefore 
maps to the ^T- a i term of Eq. Q which stabilizes the 
paramagnetic phase. In Appendix IS] we show how to 
compute the response of this term by mapping Eq. (Hi) 
to a theory of free fermions, through the subsequent ap- 
plication of the Jordan- Wigner and Bogoliubov transfor- 
mations. The resulting absorption spectrum is: 



Shop(w) = dt 



v/w 2 - (2A n )< 



9(w-2A„). (12) 



where Srii — b\bi — n. In the following we will assume 
that the Bragg spectroscopy is performed by modulating 
both the lattice and the superlattice in such a way that 
only Jfc =7r is applied. 



MI-HI phase transition 



The critical theory of the MI-HI phase transition can 
be expressed in terms of the sine-Gordon model pi) . In 
this language, the modulation of the lattice corresponds 
t<P ' ' 



r, 



hop 



St / dx cos(20) . 



(10) 



This term creates kink-anti-kink pairs, corresponding to 
particle-hole excitations. The total energy of the pair 
is given by the sum of two positive terms: the chemi- 
cal potential, A„ = 2A C , and the kinetic energy of the 
pair. As a consequence, the resulting absorption spec- 
trum S'hop(w) is non zero only for frequencies u> > 2A C . 
The relevant matrix elements can be computed using the 
form factor approach 36 39 40 , leading to the formulas sum- 
marized in Appendix [A| Close to the lowest edge of the 
continuum < u) — 2A C -C A c the spectrum follows the 
power law 



Here Q is the Heaviside theta function, and A„ = |1 — A| 
is the neutral excitation gap. As explained in SecJTT| 
the gap A„ separates the ground state from a phonon 
at momentum k = n. Lattice modulations correspond 
to a translational invariant operator and can couple only 
to excitations with zero total momentum, such as, in this 
case, a pair of counter-propagating phonons. The absorp- 
tion spectrum (12) shows a broad feature, corresponding 



to the continuum of two-particle excitations. This result 
is in contrast with the analysis of Ref.[10l which predicted 
the appearance of a sharp peak in the response to lattice 
modulations. As we will see in Sec. [V] the mean-field ap- 
proach can be fixed by a more accurate characterization 
of the elementary excitations. 

The operator corresponding to Bragg spectroscopy ftty 
is a modulation of the local operator (— l) l 5rii, which 
corresponds to the order parameter of the symmetry- 
braking phase and, in the effective Ising model, maps 
to the local magnetization af (see Eq. ([6])). In the para- 
magnetic phase (Haldane insulator), this operator cou- 
ples directly to the low-energy excitations of the phase: 
spin-flips at zero-momentum. The corresponding matrix 
element can be computed through conformal field theory 
(CFT) methodiP, leading to 



Su^iu) « |A| 2 "-M(w - A»). 



(13) 



Shop(w) ~ (u - 2A C ) 



2 A' 



(11) 



Here r\ = 1/8 is the anomalous dimension of the operator 
a x . The operator Tk =7T can additionally couple to any 
state with an odd number of spin-flips. Therefore we 



expect a continuum of three-particles excitations in Sk=n 
for w > 3|A„|. 

In the DW phase (ferromagnet) the elementary excita- 
tions are domain walls in the ferromagnetic order param- 
eter (a z ). The operator a x is local and can excite only 
an even number of domain walls. Thus, in this phase we 
expect a broad continuum, whose lowest edge is located 
at w = 2|A„|. 

Fig.s [3] and HI summarize graphically our findings for 
the absorption spectrum of lattice modulation and Bragg 
spectroscopy at the MI-HI and HI-DW phase transitions. 
The extended bright areas indicate the excitation of pairs 
or triplets of elementary excitations. The narrow absorp- 
tion peak in Bragg spectroscopy is the signature of the 
single phonon at momentum k = 7r, needed to break the 
translational invariance 



V. NON-LOCAL MEAN FIELD 

We now present an alternative method to compute 
the dynamic response of the Haldane insulator, based 
on a non-local mean field introduced by Kennedy and 
TasakP2l. This approach is valid deep in the gapped 
phase, and it complements the field-theoretical approach 
presented above, which is valid only in the vicinity of a 
phase transition. 

The first step of this method is an approximate map- 
ping between the Haldane insulator and the Haldane 
gapped phase of spin-1 chains^. The mapping can be 
made explicit by truncating the Hilbert space of the Bose 
system to three occupation states per site (for example, 
for a system with n = 1 particles per site, we keep only 
the occupation states n = 0, 1, 2 for every site). This de- 
fines an effective spin-1 model with Sf = nt — n. If we 
neglect particle-hole symmetry breaking terms, the anal- 
ogous spin-1 model is the well-known antiferromagnetic 
XXZ model 

H = E - f ( S t S r+i + K-c)+VS*S? +1 + ^(S?) 2 . (14) 



This model supports three gapped phase, the "large- U" , 
the Haldane, and the Neel phases, corresponding respec- 
tively to the MI, HI and DW phases. 

Kennedy and Tasaki 42 introduced a non-local mean- 
field theory which captures the whole phase diagram 
of the XXZ chain. Their approach relies on a unitary 
transformations^! that transforms the string operators 
into local spin operators: 
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O z string = exp ( iirJ2 S J ] S * = U- X S*U = Sf . (15) 

j<i 




FIG. 3. Energy absorption due to lattice modulations 
Sho P (oj), according to the effective field-theoretical descrip- 
tion. The horizontal axis corresponds to generic paths cross- 
ing the MI- HI (left panel) and the HI-DW (right panel) phase 
transitions. The vertical axis corresponds to the frequency 
oj. Both axes are in arbitrary units: the field-theoretical ap- 
proach describes the universal shape of the spectra around 
the phase transitions, but does not deliver quantitative infor- 
mation about their position. The value of K at the MI-HI 
phase transition S is arbitrarily set to K = 1.5. Lattice 
modulations couple to pairs of low-lying excitations near both 
transitions, leading to broad response spectra. 




FIG. 4. Energy absorption due to Bragg spectroscopy across 
the HI-DW phase transition, for the same parameters as in 
Fig. [3] Bragg spectroscopy at momentum k — -k couples to a 
single phonon whose gap closes at the HI-DW phase transition 



The operator U flips every second site with S z = ±1, 
skipping those with S z = 0. It is given explicitly bjBS 
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(16) 



The unitary operator U transforms the Hamiltonian 



( 14 1 into an unusual, but nevertheless local form 

H = UHU = -2t ]T [S*S? +1 ~ §y expfaSf + iTrS* +1 )S? +1 

i 

u 



vY / [s?sUi]+-Y,(SD 2 - 



(17) 



The transformed Hamiltonian H has an explicit Z^ X Z^ 
symmetry, generated by 7r rotations of the spins around 
the main axes. In the Haldane phase this symmetry is 
spontaneously broken and the two order parameters S x 
and S z acquire a non-zero mean value. These proper- 
ties enable a variational mean-field approximation 42 for 
the ground state. It consists in considering only site- 
factorizable wavefunctions of the form 



* 



GS . 



n i^c 



(18) 



Ki<L 



where \ip a ) =cos6»|0) ±sin9\ ±1). 

In the Haldane insulator the minimal variational en- 
ergy is obtained for cos 2 9 = (U - 2V + 4f)/(8t - 2V) 
and the ground state is four fold degenerate. As shown 
in Table U the four different states can be characterized 
by the signs of the two order parameters and reflect the 
broken Z^ X Z2 symmetry. In the DW phase cos 9 = 
and the ground state is doubly degenerate. In the MI 
phase cos 9 = 1 and the ground state is non-degenerate. 
The resulting phase boundaries, first obtained in Ref. UH 
are given in Fig. [5| 




FIG. 5. Phase diagram (solid lines) of the effective spin-1 
model (TTTJ): the three phases map to the Mott insulator (MI), 
the Haldane Insulator (HI) and the density-wave (DW) phase. 
The dashed line is one particular paths crossing the HI phase 
and used to plot Fig.s[6]and[7] 



state \ order parameter 


sign(5 a: ) 


sign(S z ) 


l^i) =cos6»|0) + sin 6»| + 1) 


+ 


+ 


|^2) =cos(9|0) -sin0| + l) 


- 


+ 


\tp 3 ) =cos0|O)+sin0|-l) 


+ 


- 


|^4) = cos0|O) -sin(9|-l) 


- 


- 



TABLE I. Fo ur m ean-field ground-states of the transformed 
Hamiltonian (171 and the respective signs of (S x ) and (S z ) 



domain wall localized at site i may be represented by the 
site factorizable wavefunction 



* 



",/3* 



Ipa] 



n 1^1 



= n 

Kj<i i+Kj'<L 

= J aa...aa pp...pp) ^ . 



(19) 



(20) 



This wavefunction can be symmetrized with respect to 
lattice translations to obtain a variational ansatz for the 
elementary excitations: 



*?*: 



E- 






* 



a,/3 1 



\aa...aa pp...pp) 



.( 21 ) 
In the Haldane insulator, the domain walls can flip either 

one of the hidden Zi symmetries or both, giving rise to 
three different excitations. Note that these excitations 
do not coincide with those predicted in Ref. [TOl As we 
will see, the present approach correctly reproduces the 
findings of the field-theoretical approach. 

A kink that flips only S x is described by the wave- 
function |*|) = |*^' 2 ) = 1 11. ..1122... 22) k (see Table 
M and is created by the untransformed operator S k . In 
terms of the original Bose-Hubbard model, this "neutral" 
excitation corresponds to a phonon (pi) with momentum 
k = 7T. In contrast, the operators ot and S^ change 
the total number of particles by one. The resulting exci- 



± 



tations | \F 
1 11...1133.r.33) 



W 



± 



n: 



* 



± 1* 



2,3 • 



j 22. ..2233. ..33 ) k correspond to a su- 



perposition of two types of kinks: 1 11. ..1133. ..33) k , 
which is a kink only in S z , and | 22. ..2233. ..33) k , which is 
a kink in both S x and S z . In the original Bose-Hubbard 
model these "charged" excitations correspond to the cre- 
ation or annihilation of one boson with momentum k. 

The energy of the three collective excitations can 
be calculated variationally: E a = (*£|ff|*g) - 
(^ GS \H\^ GS ), leading to 



A, 



El = 2*(1 + cos 2 29-2 cos 29 cos k) 



(22) 



K 



2V- 



El 
.1 — cos 2 



( 1 + cos 4 9 + 2 cos 2 9 cos k) . (23) 



We find that: (i) at the transition to the Mott insulator 
(9— >0), the energies of both excitations vanish at mo- 
mentum k — 0; (ii) at the transition to the DW phase 
(9— >7r/2), only the energy of the neutral mode (22) van- 
ishes at wavevector k = n. Both observations are in 
agreement with the numerical and field-theoretical ap- 
proaches presented above^. 



Variational wavefunctions of the excitations 



Absorption spectrum: selection rules 



The non-local mean field allows to describe the ele- 
mentary excitations of the different phases in terms of 
domain walls ( "kinks" ) of the two order parameters-* 4 } A 



Having described the elementary excitations of the 
Haldane insulator, we can now calculate the absorption 
spectra S(uj). For this task, we should compute the ma- 



trix element between the different excitations and the 
perturbation operators defined in Eq. pi and ([9| . Unfor- 
tunately, these calculations are extremely complicated, 
due to the unknown interactions between the kinks. Nev- 
ertheless, we can formulate simple selection rules that 
determine when these matrix elements can be finite and 
when they must vanish. 

To establish the selection rules, we classify all excited 
states in four topological sectors, defined by the parity of 
the number of kinks in the x and z string order param- 
eters, denoted by Af x and N z . Because the string order 
parameters of the original model correspond to the local 
spin operators of the transformed model (Eq. (fl5|)), we 
can define the Boolean numbers N x and Af z according 
to 



M a = lim sigaiSZjS^}. 



(24) 



The perturbations Mty and p)]) can connect the ground 
state only to excitations belonging to a well-defined topo- 
logical sector. In particular, the operator corresponding 
to lattice modulations, Th op , remains local after the uni- 
tary transformation ( 16 1 



UT hop TJ- 1 = Y, §jSj+i ~ Sj exp(inS* + wrSJ +1 )Sj +1 . 

(25) 
As a consequence, Thop cannot excite states with an odd 
number of domain walls, and 

( 4> a | T hop | * gs > ^ only if (AP, AP)« = (+, +) . 

(26) 
Bragg spectroscopy transforms as 

UTkU- 1 = Y,e ikxi S?l[exp(iirS?) . (27) 



j<i 



Noting that exp(MrSf) flips Sf 



-Sf, we deduce that 



T fc=7r flips N x -> -J\f x . Consequently T k=7T | * gs ) over- 
laps only with states containing an odd number of kinks 
of the x string order parameter and 

(ValZWltfg,) ^0 only if (Af x ,Af z ) a = (-,+). 

(28) 
Using these selection rules, we conclude that lattice mod- 
ulations can excite only pairs of charged or neutral ex- 
citations and therefore its absorption spectrum is domi- 
nated by two particle continua. Bragg spectroscopy Xfc =7r 
can couple to a neutral excitation at momentum k = n: 
we therefore expect a sharp resonance in the absorption 
spectrum, whose energy goes to zero at the HI-DW phase 
transition. 

The resulting absorption spectra are shown in Fig. [6] 
and [7] We find an excellent agreement with the results 
obtained from the critical theories shown in Fig. [3] and [4] 
As expected, the mean-field approach does not capture 
the precise shape of the excitation gap in the vicinity of 
the transition. It does however provide a reliable method 




FIG. 6. Energy absorption due to lattice modulations 
5ho P (w), according to the non-local mean-field description, 
along the dotted line in Fig. [5] Lattice modulations couple 
to both pairs of charged kinks (left gray area) and pairs of 
neutral kinks (right gray area). 




FIG. 7. Energy absorption due to Bragg spectroscopy, 
Sfc=7r(w), according to the non-local mean- field description, 
along the dotted line in Fig. [5] Bragg spectroscopy at k — n 
couples to single neutral kinks (white line), to pairs of charged 
kinks (gray area) 



to quantitatively predict the position of the peaks and 
continua, deep in the Haldane insulator, starting from 
the microscopic Hamiltonian M. 



VI. CONCLUSION 

In this article we described a dynamical method to 
characterize the phase transitions between three insulat- 
ing phases of interacting lattice bosons in one dimen- 
sion: the Mott insulator (MI), the Haldane insulator 
(HI), and the density- wave (DW) phases. We considered 
the response of the system to two different perturbations, 
namely lattice modulations and Bragg spectroscopy at 
momentum k = n. Measuring the energy absorption rate 
of these perturbations, it is possible to characterize the 
critical excitations of the different transitions. 
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At the MI-HI phase transition, the low energy excita- 
tions are particles and holes at zero momentum. They 
appear as a broad peak in the response to lattice modu- 
lations, and they do not couple to Bragg spectroscopy at 
k = t:. This observation indicates that the phase tran- 
sition between the Mott insulator and the Haldane in- 
sulator is of topological origin: the critical excitations 
at a topological phase transition trivially commute with 
all the symmetries of the problem, i. e. translational 
invariance and particle conservation. 

At the HI-DW phase transition, the low energy exci- 
tations are phonons at the edge of the Brillouin zone. 
These excitations appear as a narrow peak in the re- 
sponse to Bragg spectroscopy at k — tt. They hold a 
non-trivial quantum number (momentum) , showing that 
the HI-DW phase transition is not topological, and spon- 
taneously breaks a local symmetry (translational invari- 
ance). 

These qualitative features are confirmed by the field- 
theoretical analysis of the phase transitions. More quan- 
titative predictions are obtained using a non-local mean- 
field approach. In contrast to the analysis of Ref. [TQl the 
present approach correctly reproduces the shape of the 
spectra predicted by the field-theoretical calculations, at 
both the MI-HI and HI-DW phase transitions. The non- 
local mean-field can be therefore used to interpolate be- 
tween the two phase transitions, and is expected to pro- 
vide reliable results deep in the topological phase (HI), 
where the field theories cease to apply. This method does 
not rely on effective parameters of the continuum theory, 
which are in general not known, but rather applies di- 
rectly to the microscopic Bose-Hubbard model. 

The techniques developed here can be extended to 
other types of dynamic probes, such as Bragg spec- 
troscopy at a general momentum k, local perturbations, 
small quenches^, as well as to quantum phase transitions 
of other systems. Among those, we would like to briefly 
discuss quantum systems of interacting fermions in one 
dimension. The extended Bose-Hubbard model consid- 
ered here can be approximately mapped into a fermionic 
problem, by splitting the single chain at integer filling 
into two chains at half-integer filling, and then applying 
the Jordan- Wigner transformation. This transformation 
maps the local occupation rij and the nearest neighbor 
tunneling b\bi + i of the bosons into the corresponding op- 
erators of the fermions. As a consequence, we expect the 
absorption spectrum of fermions to be identical to the 
one of bosons. The topological phase transition between 
the Mott insulator and the Haldane insulator maps into a 
phase transition between two topological phases of inter- 
acting fermions in one dimension, belonging to the more 
general Z s group of Ref. 47. (The fermionic model pos- 
sesses additional symmetries which map into local sym- 
metries of the bosonic model.) It would be interesting 
to extend this result to the transitions to the other six 
phases. 

Field-theoretical techniques should allow to extend the 
present approach to topological phases of fermions and 



bosons in higher dimensions. Recently, field theories de- 
scribing gapped phases of fermions^ and bosons^ have 
been found. However, less is known about the critical 
theory of the quantum phase transitions between the dif- 
ferent topological phases. In the case of topological band 
insulators, it is natural to expect these critical theories to 
be non-interacting (Gaussian), in analogy to the present 
case. In general it should be possible to compute the re- 
sponse to physical probes by expressing the correspond- 
ing perturbations in terms of the critical theories. One 
should than be able to experimentally probe that the crit- 
ical excitations of the topological phase transition respect 
all the symmetries of the systems. 



VII. ACKNOWLEDGMENT 

We wish to thank E. Berg, T. Giamarchi, S. Huber, S. 
Kivelson for many useful discussions. Part of the results 
presented in this work has been submitted to the Fein- 
berg Graduate School, as a chapter of the Ph.D. thesis 
of the author, under the supervision of E. Altman. 



Appendix A: Form-factor approach 



According to the form-factor approac h 1 36 ! 39 ! 40 !, the 
spectrum of lattice modulations is given by: 



S'hop(w) 



F cos [6(u,g)} 



where A c is the charge gap ("soliton mass") and 



(Al) 




2A ' V 4A : - ; 
cosh(<9/2) 



- I 



sinh 



5^»-«rJ 



F(0), 



F(9) = sinh(0/2) exp / dx K{6,x) , 



2 . 



K{6,x) 



(sin [x{6 + iTt)/2]) sinh [(tt - j)x/2] 
x sinh(7rx) cosh(7ra;/2) sinh(7x/2) 
itK 



2-K 



(A2) 
(A3) 

(A4) 

(A5) 
(A6) 



Appendix B: Jordan- Wigner and Bogoliubov 
transformations 



In this appendix we compute the absorption spectrum 
of Thop = V\ S" through the subsequent application of 
Jordan- Wigner and Bogoliubov transformation. Apply- 
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ing Jordan- Wigner transformation to p} we obtain: 

Rising = Y, ^ ~ X ^i + ^) (4+1 + 4>i+l) (Bl) 

i 

= y^(l - Acos(fc))V'fcV'fe + Asin(fc)^fcV'fe + H.c. 

Here u = 2Aa and A = 2(A — 1). We now diagonalize 
( |B1[ ) by the Bogoliubov transform: 7^. = cos(0fc/2)V>& — 
ism(9k/2)tjj_ k , where tan(0fe) = k/A. 

We now express a z in terms of Bogoliubov quasi- 
particles creation and annihilation: 

J2 < = sin(0 fe ) 7 i7i fe + cos(» fc )7hfc + Hc - ( B2 ) 

i 

The required matrix element is: 

M M = ( Vfc I 5Z ^ I *gs ) = sin (^) ~ — • ( B3 ) 



Here \ipk) is a state with two Bogoliubov quasi-particles 
(at momentum k and —k) and energy 2gfc = 2\/k 2 + A 2 . 
The computation of the absorption spectrum is straight- 
forward: 



ShopM « ^ |M M | 2 <5(w - 2e fe ) 



0fc 



* / w=2e fc 



ft 



V^ 2 - (2A) 2 



9(w-2A) (B4) 



Here <d(u) is the Heaviside step function. 
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